program mynonlin_lf
	version 8.1
	args lnf theta1 theta2 theta3 sigma
	quietly replace `lnf' = ///
	ln(normden($ML_y1,`theta1'+`theta2'*X3^`theta3',`sigma'))
end
